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Abstract 

We investigate the rank of the adjacency matrix of large diluted random graphs: for 
a sequence of graphs (G„)„> converging locally to a Galton- Watson tree T (GWT), we 
provide an explicit formula for the asymptotic multiplicity of the eigenvalue in terms of the 
degree generating function ip* of T. In the first part, we show that the adjacency operator 
associated with T is always self-adjoint ; we analyse the associated spectral measure at 
the root and characterize the distribution of its atomic mass at 0. In the second part, we 
establish a sufficient condition on <y9* for the expectation of this atomic mass to be precisely 
the normalized limit of the dimension of the kernel of the adjacency matrices of (G„) n >rj. 
Our proofs borrow ideas from analysis of algorithms, functional analysis, random matrix 
theory, and statistical physics. 
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1 Introduction 

In this paper we investigate asymptotical spectral properties of the adjacency matrix of large 
random graphs. To motivate our work, let us briefly mention its implications in the special case 
of Erdos-Renyi random graphs. Let G n = (V n ,E n ) be an Erdos-Renyi graph with connectivity 
c > on the vertex set V n = {1, • • • , n}. In other words, we let each pair of distinct vertices ij 
belong to the edge-set E n with probability c/n, independently of the other pairs. The adjacency 
matrix A n of G n is the n x n symmetric matrix defined by (A n )ij = l((ij) € E n ). Let \\{A n ) > 
■ • • > A ra (vl n ) denote the eigenvalues of A n (with multiplicities) and 

1 - 

i=l 

denote the spectral measure of A n . Our main concern will be the rank of A n \ 

rank(^4„) = n — dimker(A n ) = n — nfi n ({0}). 



Theorem 1 (i) There exists a deterministic symmetric measure fx such that, almost surely, 
for the weak convergence of probability measures, 

lim [i n = /x. 

n— >oo 

(ii) Let < q < 1 be the smallest solution to q = exp(— cexp(— cq)). Then almost surely, 

lim Mn ({0}) = /i({0}) = q + e- c « + cqe~^ - 1. 
n— >oo 

In other words, almost surely, 

rank(^4„) „„ 
lim = 2-q- e~ cq - cqe~ cq . 1 

n— too fi 
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Apart from an improvement of the convergence, part (i) is not new : the convergence in prob- 
ability was first rigorously proved by Khorunzhy, Shcherbina and Vengerovsky [16] ; for an 
alternative proof, see [TU]. Note that it only implies limsup n /i n ({0}) < fi({0}). 

In the sparse case, i.e. when the connectivity c grows with n like alogn, the rank of A n has 
been studied by Costello, Tao and Vu [13] and Costello and Vu |12j . Their results imply that 
for a > 1, with high probability dimker(A n ) = while for < a < 1, dimker(j4 n ) is of order of 
magnitude n 1_cl . Our theorem answers one of their open questions in |12| . 

The formula ([TJ already appeared in a remarkable paper by Karp and Sipser [15] as the 
asymptotic size of the number of vertices left unmatched by a maximum matching of G n . To 
be more precise, the function G i— > dimker(G) is easily checked to be invariant under 'leaf 
removal', i.e. if G' is the graph obtained from G by deleting a leaf and its unique neighbor, 
then dimker(G') = dimker(G). Karp and Sipser [15] studied the effect of iterating this leaf 
removal on the random graph G n until only isolated vertices and a 'core' with minimum degree 
at least 2 remain. They showed that the asymptotic number of isolated vertices is approximately 
(2 — q — e~ cq — cqe~ cq )n as n — > oo, and that the size of the core is o(n) when c < e. Thus, 
([U follows by additivity of G H> dimker(G) on disjoint components, as observed by Bauer and 
Golinelli [6]. However for c > e, the size of the core is not negligible and the same argument 
only leads to the following inequality : 

,. . „dimker(A„) _„„ 

hm mf >q + e cq + cqe cq - 1. 

n— >oo n 

Bauer and Golinelli [6] conjectured that this lower bound should be the actual limit for all c, 
which is equivalent to saying that asymptotically the dimension of the kernel of the core is zero. 
The proof of this conjecture follows from our work (see Section 

Our results are not restricted to Erdos-Renyi graphs. They will in fact hold for any sequence 
(G n ) n >i of random graphs converging locally to a rooted Galton- Watson tree (GWT), provided 
the latter satisfies a certain degree condition. The precise definition of local convergence is 
recalled in Section [3l It was introduced by Benjamini and Schramm [7], Aldous and Steele 
[3]. A rooted GWT (see [2]) is characterized by its degree distribution i 7 *, which can be any 
probability measure with finite mean on N : the root has offspring distribution F* and all 
other genitors have offspring distribution F, where for all k > 1, F(k — 1) = kF^i^k)/ ^2 i £F if (£). 
In the case of Erdos-Renyi graphs with connectivity c, the limiting tree is simply a GWT with 
degree distribution F* = Poisson(c). 

The adjacency operator A of a GWT T = (V,E) is a densely defined symmetric linear 
operator on the Hilbert space l 2 (V) defined for i,j in V by 

(Aei,ej) = 1 (ij € E) , 

where for any i € V, e\ denotes the base function j £ V h- > l(j = i). As we will show, if 
has a finite second moment, then A has almost surely a unique self-adjoint extension, which 
we also denote by A. Consequently, for any unitary vector ip 6 Dom(^4), the spectral theorem 
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guarantees the existence and uniqueness of a probability measure fi^p on R, called the spectral 
measure associated with such that for any k > 0, 

(A k il>,ifj) = / x k dn i3 {x). 

In particular, we may consider the spectral measure /ir associated with the vector e@, where 
is the root of the rooted tree T. Our first main result is an explicit formula for E/x*r({0}), the 
expected mass at zero of the spectral measure at the root of a rooted GWT T. 

Theorem 2 Let T be a GWT whose degree distribution has a finite second moment, and let 
tp* be the generating function of F*. Then, E//t({0}) =max 3 . e rQ jl iM(x), where 

M(x) = <p^(l)xx + (p*(l-x) + <p*(l — x)-l, with T = (p , ;¥ (l — x)/<d ¥ (l). 

In the special case of regular trees, the measure \it can be explicitely computed and turns 
out to be absolutely continuous, so /iT({0}) = 0. In contrast, one may construct GWT's with 
arbitrary large minimum degree and such that E/xt({0}) > 0. The following example is taken 
form [9 J and is due to Picollelli and Molloy: set d > 3 and take cp*(x) = j^x d + j^x d3 . Figure 
[H gives a plot of M for the case d = 3, showing that E/u^({0}) > in this case. 




When is a Poisson distribution with mean c, the corresponding quantity max^grou M(x) is 
precisely (pQ), and it already appeared in Zdeborova and Mezard [HI Equation (38)], as a 'cavity 
method' prediction for the limiting fraction of unmatched vertices in a maximum matching. 

To the best of our knowledge, the formula was unkwown for general GWT's. However, Bauer 
and Gollineli |5j have computed explicitly the asymptotic rank of the uniform spanning tree on 
the complete graph of size n. Also Bhamidi, Evans and Sen [8] have recently analyzed the 
convergence of the spectrum of the adjacency matrix of growing random trees. 

Our second main result (Theorem [13]) states that for any sequence of random graphs (G n ) n >o 
converging locally in distribution to a GWT, we have lim n n _1 rank(^4 n ) = 1— E//t({0}), provided 
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the first local extremum of the above function x i— > M(x) is a global maximum on [0,1]. We 
have left open the case where the global maximum of M is not the first local maximum (see 
Section 0}. 

Our detailed analysis of the atomic mass at of the limiting spectral measure \x remains only 
a small achievement for the global understanding of this measure. For example, for Erdos-Renyi 
graphs, the atomic part of [i is dense in K and nothing is known on the mass of atoms apart 
0. There is also a conjecture obout the absolutely continuous part [i ac of the measure \x : we 
say that \x has extended states (resp. no extended state) at E 6 M if the partition function 
x i — y [i ac (—oo,x) is differentiable at x = E and its derivative is positive (resp. null). This 
notion was introduced in mathematical physics in the context of spectra of random Schrodinger 
operators ; a recent treatment can be found in Aizenman, Sims and Warzel [I]. For Erdos-Renyi 
graphs, Bauer and Gollineli have conjectured that \x has no extended state at E = when 
< c < e, and has extended states at E = when c > e. More generally, one may wonder 
whether [i ac = when < c < e. Finally, the existence of a singular continuous part in \i is 
apparently unkown. 

The remainder of the paper is organized as follows: in Section [2 we analyze the adjacency 
operator of a GWT. In Section [2T3l we study ht({0}) and prove Theorem [2j In Section El we 
prove finally the convergence of the spectrum of finite graphs and the convergence of the rank. 
The proof of Theorem [1] is given in Appendix [5j 

2 Locally finite graphs and their adjacency operators 

A rooted graph is the pair formed by a graph G with a distinguished vertex € V, called the 
root. There is a canonical way to define a distance on V: for each u, v € V, the (graph)-distance 
is the minimal length of a path from u to v, if any, and oo otherwise. For a rooted graph G with 
root and t an integer, we will denote by (G)t the rooted subgraph spanned by the vertices at 
distance at most t from the root. In all this section, we consider a locally finite rooted graph 
G = (V, E) with root denoted by 0. 

2.1 Adjacency operator 

Consider the Hilbert space 



Denote by Hq C £ 2 (V) the dense subspace of finitely supported functions, and by (ej) igl/ the 
canonical orthonormal basis of £ 2 (V), i.e. e\ is the coordinate function j € V i— > l(i = j). By 
definition, the adjacency operator A of G is the densely-defined linear operator over £ 2 (V) whose 
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domain is Hq and whose action on the basis vector e\, i G V is: 

Ae i = Y e i- 

Note that Ae K G £ 2 (F) since G is locally finite. Moreover, for all i,j G 

(Aei.ej) = l{ij G -B} = (Aej.ei). 

Therefore, the operator y4 is symmetric, and we may now ask about the self-adjointness of its 
closure, which is again denoted by A. The answer of course depends upon G, but here is a 
simple sufficient condition that should suits all our needs in the present paper. 

We define the boundary of a subset 5 C V as dS = {ij € E: i G S,j ^ S}, and the boundary 
degree A(dS) as the maximum number of boundary edges that are adjacent to the same vertex. 



Proposition 3 For A to be self-adjoint, it is enough that V admits an exhausting sequence of 
finite subsets with bounded boundary degree: 

(A) There exists finite subsets Si, S2, ... QV such that U n S n = V and sup n A(dS n ) < 00. 

Proof. Denote by A* the adjoint of A. By the basic criterion for self-adjointness (see e.g. Reed 
and Simon, Theorem VIII. 3, |18j), it is enough to show that is the only vector ip G Dom(A*) 
satisfying A*ip = ±iip. Consider such a tjj (let us treat, say, the +i case), and define the following 
flow along the oriented edges of G: 



(i->j) = 3(^(i)V(j)J =-(j->i), 
for all ij G E. The amount of flow created at vertex i G V is then: 

Y (i->j) = 3U(i) Y, OT] =^W(i)ei),V} = 9^(i)ei,^» = |^(i)| 2 . 

Now, by anti-symmetry of the flow, the total amount of flow created inside any finite subset 
5cy must equal the total amount of flow escaping through the boundary OS: 

Ym)\ 2 = E(^j)- 

i€S ijedS 

Therefore, using (i — > j) < |^(i)||V ; (j)| and twice the Cauchy-Schwartz inequality, we find : 

£hKi)i 2 < ( E i^)i 2 E ( E M)\ 

ies \ieds- ieds- \jeS c rWi 



_> 



< A(dS) [ Y I^WI 2 E l^(j)' 2 
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where we have written N{ for the set of neighbors of i, dS~ and dS + for the sets of vertices 
dS n S and dS n S c respectively. Finally, take S = S n , and let n —¥ oo: the exhaustivity 
Un'S'n = V ensures that the left-hand side tends to YlieV I^WI 2 = IIV'II 2 an d also that : 

£ |V>(i)| 2 --^0and \m\ 2 —^0- 
Since sup n A(dS n ) < oo, the right-hand side vanishes, and we obtain the desired ||V>|| = 0. □ 
2.2 Spectral measure 

We now assume that the adjacency operator A is self-adjoint. The spectral theorem then 
guarantees the validity of the Borel functional calculus on A : any measurable function / : R — >■ 
C may now be rigorously applied to the operator A just as one would do with polynomials. 
Denoting by the spectral measure associated with the vector e , we may thus write : 

(f(A)e ,e ) = [ f(xW G (x) (2) 

JR 

for any / G £c(mg)- Taking f(x) = x n (n <G N), we obtain in particular 

7„ = (A n e , e ) = J x n dfiG(x) = # {paths of length n from to in G}. (3) 

Since 1 1 e<2> 1 1 = 1> the spectral measure [It is a probability measure on R. We will now study 
its Cauchy-Stieltjes transform. By definition, the Cauchy-Stieltjes transform of a probability 
measure (i on R is the holomorphic function defined on the upper complex half-plane C + by 

f dp(x) 

inn,: zv^ I . 

Jr x - z 

Note that m M belongs to the set % of holomorphic functions / on C + satisfying 

Vz € C+,3/(z) > and \ f(z)\ < (^z)-\ 

which is compact in the normed space of holomorphic functions on C + (Montel's Theorem). 

Henceforth, we will assume that G is a rooted tree T. We write j y i to mean that i G V is 
an ancestor of j € V, and we let Tj be the subtree of T restricted to {j G V,j >z i}, rooted at i. 
Its adjacency operator A\ is the projection of A on Vect (ej, j >z i). Since it is also self-adjoint, 
we may consider its spectral measure associated with the vector e- u and its Cauchy-Stieltjes 
transform my.. The recursive structure of trees implies a simple well-known recursion for the 
family {m T .) ieV : 
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Proposition 4 The family (mT^icv is solution in T~L V to the system of equations, for all z 6 
C+, 

M*) = -\z+ E /jto] . ( 4 ) 

\ je-D(i) / 

where D(i) = {j >- i, [j| = |i| + 1} denotes the set of immediate children of i. 

Proof. As we will see, the recursion follows from a classical operator version of Schur comple- 
ment formula (see e.g. Proposition 2.1 in Klein [T7] for a similar argument). We write the proof 
for completeness. Define the operator U on £ 2 (V) by its matrix elements: 

(£7e ,ei) = (Ue h e ) = 1 = (Ae ,ei), 

for all i € D{&), and (C/ej^k) = otherwise. We then have the following decomposition: 

A = U+ A i9 

ie£>(0) 

where Ai, is the projection of 4 on V; = Vect(ej,j >: i). Since A and A = ©i € D(0) are 
self-adjoint operators, their respective resolvents 

R: z ^ (A - ziy 1 , R: z ^ (I- zl)~ l 

are well-defined on C+, and the resolvent identity gives 

R{z)UR(z) = R(z) - R{z). (5) 

In particular, for all k € V, 

(R(z)UR(z)e ,e ]i } = (R(z)e , e k ) - {R(z)e , e k ). 
Now, using the definition of U, we may expand the left-hand side as 

(R(z)e ,e k ) ^2 i R ( z ) e 0i e i) I + ( R ( z ) e 0, e 0) E (^( z ) e i' e k) I ■ 
ieD(0) J \ ieD(0) J 

But R(z)e = —z~ l e and each Vj,i € -D(0), is stable for R. Therefore, in the special case 
where k = 0, the above equality simplifies into : 



- E ( R ( z ) e 0' e O = ( R ( z ) e 0> e 0) + -> 



Z ieD(0) 



while for k £ D(0), it gives : 

(R(z)e ,e }(R(z)e k ,ek) = (R(z)e ,e k ). 
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Combining both, we finally obtain 

(R(z)e ,e ) = - iz+ ^ (R{z)e h e { ) 

\ ie-D(0) 

which, by (|2|) with /(x) = (x — z) -1 , is precisely 

m T {z) =-\z+ Yj m Ti{ Z ) 
V ieD{0) 



□ 



When T is finite, the set of equations ([!]) uniquely determines the Cauchy-Stieltjes transforms 
(myJigV) which can be computed iteratively from the leaves up to the root. Under an extra 
conditon on T, this extends to the infinite case. Recall that (T) n denote the truncation of T to 
the first n generations. In what follows, we will make the additional assumption: 

(B) limsup |<9(T) n | « < oo. 
n— too 

Proposition 5 If T satisfies Assumption (B), then (mrJieV is the unique solution in H v to 
the system of equations §4ty, and for all i £ V , 

in the sense of compact convergence on C_|_. 

Proof. If (/i)i g y € H and (<?i)j e y € H v are solutions to the system of equations (jl]), then 
we can write, for all i G V, z € C+, 



\Si{z)-9i{z) 



£j€r>(i)(/j(z) ~9j(z)) 



v v " jeD(i) 



v * + EjGD(i) /j (*) J + EjeD(i) sa (*) y 

Iterating this n times, and then using the uniform bound \fj(z) — gj(z)\ < 2 we obtain 

- *(*)\ < E l/j(')-^)l< ( S - 

v v '' ied(Ti) n K v >> 

Therefore, we see that under assumption (B), 

Vie VAfi(z) - 9i (z)\ =0 

as soon as is sufficiently large, hence for all z € C+ by holomorphy. Finally, denote by 
M n the denumerable vector of holomorphic functions (wi(in n )iev S Since % is compact, 
the sequence (M ra ) n >o is relatively compact, and since each vector M n satisfies the partial set 
of equations (jl]) corresponding to the truncated tree (T) n , any limit point must satisfy the 
global set of equations (jl]) corresponding to the full tree T, so Mqo is nothing but (mrjiey- 
Therefore, the sequence of vectors (M n ) n >o converges to M m = (myjigy, and this is exactly 
(El). ' □ 
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2.3 Atomic mass at zero 



Our goal here is to characterize ht({0}), the atomic mass at zero of the spectral measure ht- 

Proposition 6 If T satisfies Assumption (B), then, the family (//r.({0}))i e y is the largest so- 
lution in [0,1]^ to the system of equations 

!+ E ( E **) ) > ( 7 ) 

\keD(j) / / 

with the conventions 1/0 = oo and l/oo = 0. 

Proof. Since T is acyclic, ([3]) ensures that the measures ^T i; i G V are symmetric. Therefore, 
for alU > 0, i G V 

m Ti (it) = / k ^d m (x) + i / k -r^M*) = i / k ^^t ; (x). 
Hence, if we define hr^t) := —it my. (it) € [0, 1], then by the dominated convergence theorem, 

M«> = - -,({0}). 

But, iterating once equation we get 

hr i (t) = (l+ E (* 2 + E M*)] ] > ( § ) 

\ jeD(l) V kGD(j) / / 

so that letting t — > yields exactly that (/^.({OD^gv" must satisfy ([7]). 

Again, when the rooted tree T is finite, this recursion characterizes the family (MT ; ({0})) ie y, 
since it can be computed iteratively from the leaves up to the root. However, when T is infinite, 
([7]) may admit several other solutions. Fortunately, among all of them, (/^T;({0})) ig y is always 
the largest. To see why, consider any solution (xi)i^T £ [0j 1]^- Fixing t > 0, let us show by 
induction that for all n G N, 

Vi G V, Xi < h {Tl ) 2n (t) := -itm {T . hn (it), (9) 

This will conclude our proof since we may then let n — > oo to obtain x\ < (t) by Proposition 
[5j and let finally t — > to reach the desired x\ < ^x-({0}). The base case n = is trivial because 
the right-hand equals 1. Now, if ([9]) holds for some n G N, then for all i G V, 

xi = ( 1 + E ( E Xk 

< (i+ E h 2 + E h (T kh jt)\ ) =h (Tihn+2 (t), 
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where the first equality follows from the fact that (xi)i^x satisfies ([7]), the middle inequality 
from the induction hypothesis, and the last equality from (|SJ) applied to {T{)2n+2- 



□ 



2.4 Galton- Watson trees 

We now apply the above results to Galton- Watson trees. Let be a distribution on N with 
finite mean, and let T be a GWT with degree distribution F*, i.e. a random locally finite rooted 
tree obtained by a Galton- Watson branching process where the root has offspring distribution 
F* and all other genitors have offspring distribution F, where 



In the rest of this paper, we will make the following second moment assumption on the distri- 
bution F*: 52 k k 2 F*(k) < oo, or equivalently Ylk kF(k) < oo. It is in fact a sufficient condition 
for all the previous results to hold almost surely : 

Proposition 7 If F* has a finite second moment, then T satisfies (A) and (B) with probability 
one. In particular, the adjacency operator A is almost surely self-adjoint, and the atomic mass 
at zero of the spectral measure at the root of T is characterized by the fixed-point equation (71). 

Proof. Let N denote a generic random variables with law F. For (B), it is well-known (and 
easy to check by a martingale argument) that the size of the n th generation in a GWT with 
offspring distribution F behaves like E n iV as n — >• oo, in the precise sense that almost surely, 
n -1 log |<9(T) n | — > EiV, which is finite by assumption. As far as (A) is concerned now, if T is 
finite there is nothing to do. Now if T is infinite, we build an exhausting sequence of finite vertex 
subsets with uniformly bounded boundary degree as follows: the finite first moment assumption 
on F guarantees the existence of a large enough integer k > 1 so that 



for each vertex of T, color it in red if it has less than k children and in blue otherwise. If the 
root is red, set Si = {0}. Otherwise, the connected blue component containing the root is 
a GWT with average offspring J2k>K kF(k) < 1, so it is almost-surely finite and we define S± 
as the set of its vertices, together with their (red) external boundary vertices. Now for each 
external boundary vertex i € dSf, we repeat the procedure on the subtree T\, and we define 
52 as the union of S\ and all the resulting subsets. Iterating this procedure, we obtain an 
exhaustive sequence of subsets S\ , S2 , ■ ■ ■ C V whose boundary degree satisfies by construction 
A(dS n ) = k, which is exactly (A). □ 

Owing to the recursive distributional nature of GWTs, the set of equations ([7]) defining 
Mt({0}) takes the much nicer form of a Recursive Distributional Equation (RDE), which we 




(10) 



e 




(11) 
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now explicitate. We denote P(N) (resp. FQO, 1])) the space of probability distributions on N 
([0,1] resp.). Given F,F' G V(N) and v G V ([0, 1]), we denote by Q F ,F'(v) the distribution of 
the [0,1]— valued r.v. 

y = — / ... — — (12) 



where N ~ F, N/ ~ i 7 ', and JQj ~ z/, all of them being independent. With this notation in 
hands, the previous result implies the following : if F* has a finite second moment, then /^({O}) 
has distribution Of*,f(^o)' where F is given by (I10p and v$ is the largest solution to the RDE 

"o = 8vK)- (13) 



The remaining of this section is dedicated to solving (113j) when F* has a finite second moment. 
We will assume that F*(0) + F*(1) < 1, otherwise F = 5o and z/q = 5i is clearly the only solution 
to ([13]). We let tp*(z) = J2 n >o F*{ n ) zTl be the generating function of F*. For any x G [0, 1], we 
set x = </?*(! — X )/V*(l) an d we define 

M(x) = <pi(l)xx + (p*(l - x) + </>*(l - x) - 1. 

Observe that M'(x) = </?"(l — x) (x — x) and therefore, any x G [0, 1] where M admits a local 
extremum must satisfy x = x. We will say that M admits an historical record at x if x = x 
and M(x) > M(y) for any < y < x. Since [0, 1] is compact and M is analytic, there are only 
finitely many such records. In fact, they are in one-to-one correspondence with the solutions to 
the RDE (HHP. 



Theorem 8 If pi < ... < p r are the locations of the historical records of M , then the RDE 
U3\) admits exactly r solutions ; moreover, these solutions can be stochastically ordered, say 
v\< . . . <v r , and for any i G {1, . . . , r}, 

(%) u l ({o} c )=Pi ; 

(ii) Si^iK^i) has mean M(pi). 
In particular, E [/^r({0})] = maxa;e[0,i] M(x). 

It now remains to prove Theorem The space V ([0, 1]) is naturally equipped with 

- a natural topology, which is that of weak convergence: 

fj, n > \i ItpdfAn > Itpdn for any continuous function ip: [0, 1] — > R; 

n— ¥oo J n— ¥00 J 

- a natural order, which is that of stochastic domination: 

Ml < M2 Jfdfii < Jfdfi2 for any continous, increasing function ip: [0, 1] — > R. 
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The proof is based on two lemmas, the first one being straightforward. 



Lemma 9 For any F,F' G 7 ? (N)\{5o}, @f,f' i> s continuous and strictly increasing on V ([0, 1]). 



Lemma 10 For any v G V ([0, 1]), letting p = v ({0} c ), we have 

(i) Q FtF (u)({0} c )=f 

(ii) if ®f,f{v) < v, then the mean of ®F*,F{y) is at least M(p). 
(Hi) if ®f,f{ v ) > v, then the mean of &f*,f(v) is at most M(p); 

In particular, if v is a fixed point of &f,f, then p = p and @f,,f(v) has mean M(p). 

Proof. In equation (|12p it is clear that Y > if and only if for any i G {1, . . . , N}, there exists 
j G {1, . . . , N/} such that Xy > 0. Denoting by ip the generating function of F, this rewrites: 

eF,F(^)({o} c ) = ^(i-^(i-^({o} c ))). 

But from (fTUI) it follows that <£>(•) = </?*(-)/(^(l), i.e. (p (1 — x) = W, hence the first result. 

Now let X ~ v, N* ~ F*, N ~ F, and let S,Si,... have the distribution of the sum of N 
i.i.d. copies of X, all these variables being independent. Then, Qf„,f( 1 ') has mean 



E 



1 



E 



iV» c -1 



^(l-p)-^(l)E 
p*(l- (1)E 



— ■ J-{Vj=l...Af*,S' 4 >0} 
5- 1 



s- 1 + 1 + Ef=i « 
y 

yTs 1{5>0} 



JV c.-l {5>0,Vj=l...Af*,5 4 >0} 



where the second and last lines follow from (|10p and Y ~ O^?^^), respectively. Now, for any 
s > 0, x i-> is increasing and hence, depending on whether Qf,f{^) > v or < ^, 

®f*,f{v) has mean at most/least: 



^(l-p)-^(l)E 



L {5>0} 



^(l-p)-M(l)E 



1 



1 + iV 



-{7V>1} 



A' 



with jV = ^l {Xi>0} . 

i=l 

(14) 

But using the definition (jlOp and the well known identity (n + 1) = (<i + 1) (^i ) , one can 
easily check that : 



1 

1 + iV 



L {iV>l} 



p.(i-p)-M(i)E f ( n )E 



n>l 



d=l 



d d+1 



M{p). 
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□ 

We now have all the ingredients we need to prove Theorem 

Proof. Let p £ [0, 1] such that p = p, and define vq = Bernoulli (p). iFrom lemma [TU1 we know 
that @f,f (^o) ({0} c ) = Pi an d since Bernoulli(p) is the largest element of "P([0, 1]) putting mass 
p on {0} c , we have @ff (^o) < ^o- Immediately, Lemma [9] guarantees that the limit 

h>oo = lim \ Q k FF (z^ ) 

k— >oo 

exists in 7 ? ([0, 1]) and is a fixed point of Qf,f- Moreover, by Fatou's lemma, the number 
Poo = ^oo ({0} c ) must satisfy p^ < p. But then the mean of Of*,f(^oo) must be both 

- equal to M(poo) by lemma [TU1 with v = ; 

- at least M(p) since V/c > 0, the mean of Of*,f (©ff^o)) is at least M(p) (lLemma [TU1 
with v = @ k FF (u )). 

We have just shown both M(p) < M^poo) and p^ < p. From this, we will now deduce the 
one-to-one correspondence between historical records of M and fixed points of Gf,f- We treat 
each inclusion separately: 

- If M admits an historical record at p, then clearly p^ = p, so is a fixed point satisfying 

^oo({0} C )=F- 

- Conversely, considering a fixed point v with v ({0} c ) = p, we want to deduce that M 
admits an historical record at p. We first claim that v is the above defined limit v^. 
Indeed, v < Bernoulli(p) implies v < (0f,f is increasing), and in particular p < p^. 
Therefore, p = p^ and M(p) = M(p OQ ). In other words, the two ordered distributions 
©f*,f( i/ ) < ©f*,f (^oo) share the same mean, hence are equal. This ensures v = v^. Now, 
if q < p is any historical record location, we know from part 1 that 

Aoo = lim \ @ FF (Bernoulli (g)) 

fc— >oo ' 

is a fixed point of 0f,f satisfying ({0} c ) = q. But q < p, so Bernoulli^) < Bernoulli^), 
hence < v^. Moreover, this limit inequality is strict because A^ ({0} c ) = q < p = 
^oo({0} c ). Consequently, 0f*,f(Aoo) < ©f^f^oo) and taking expectations, M{q) < 
M(p). Thus, M admits an historical record at p. 

□ 
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3 Convergence of the spectral measure 



3.1 Local convergence of rooted graphs 

In this paragraph, we briefly recall the framework of local convergence introduced by Benjamini 
and Schramm [7], Aldous and Steele [3], see also Aldous and Lyons [2]. 

We recall that for integer t, (G)t is the rooted subgraph spanned by the vertices at distance 
at most t from the root. We consider the set G* of all locally finite, connected rooted graphs, 
taken up to root-preserving isomorphism. With the terminology of combinatorics, G* is the set 
of rooted unlabeled connected locally finite graphs. We define a metric on G* by letting the 
distance between two rooted graphs G% and G2 be 1/(1 + T), where T is the supremum of those 
t > such that there exists a root-preserving isomorphism from (G\)t to (G^)*- Convergence 
with respect to this metric is called local convergence. 

This makes G* into a separable and complete metric space, see §2 in [2]. In particular, we 
can endow G* with its Borel cr-algebra and speak about weak convergence of random elements 
in G* ■ Specifically, a sequence of probability distributions pi , P2 , . . . on G* converges weakly to 
a probability distribution p, denoted by p n => p, if 



for all bounded continuous function / : G* — > K. This is called the local weak convergence. 

Let us finally mention three important examples of random graph sequences that converge 
locally weakly to Galton- Watson trees. The Erdos-Renyi graphs with connectivity c on the vertex 
set {1, . . . , n}, rooted at = 1 converges locally weakly to the GWT with degree distribution 
Poisson(c). The uniform /c-regular (k > 2) graph on {l,...,n}, rooted at = 1, converges 
weakly to the infinite /c-regular tree. More generally, if is a degree distribution on N with finite 
mean, the random graph-sequence with asymptotic degree distribution F* converges to the GWT 
with degree distribution F*. Note that in the above examples, the vertices are exchangeable and 
the choice = 1 is arbitrary: equivalently, we could have chosen uniformly at random among 
all vertices, independently of the edge structure. 

3.2 Continuity of the spectral measure 

Since the elements of Q* have countably many vertices and are only considered up to isomor- 
phism, we may without loss of generalities embed all vertices into the same, fixed generic vertex 
set V, say the set of finite words over integers: the root is represented by the empty- word 0, 
and vertices at distance t from the root are represented by word of length t in the usual way. All 
adjacency operators can thus be viewed as acting on the same Hilbert space £ 2 (V), their action 
being defined as zero on the orthogonal complement of the subspace spanned by their vertices. 
Note that this does not affect the spectral measure at the root px- 
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If (G n ) is a converging sequence in G*, say to G £ G*, we may even relabel the vertices in 
a consistent way so that the root-preserving isomorphisms appearing in the definition of local 
convergence become identities: for every t £ N, there exists nt £ N such that 

n>n 4 (G n ) t = (G) t . (15) 

Fixing a word i £ V, and setting t equal 1 plus the distance from i to the root above, we 
obtain that for all n > nt, i is a vertex of G n if and only if it is a vertex of G, and in that case 
its neighbors in G n are exactly its neighbors in G. In other words, A n e\ = Ae\. By linearity, it 
follows that any finitely supported vector i/j : V — > C must satisfy 

e 2 (v) 
A n ifj — ^> A0, 

n— >oo 

and since those ip are dense in £ 2 (y), Theorem VIII. 25(a) in Reed and Simon [18] guarantees 
that A n — > A in the strong resolvent sense, provided of course that A,A\,... are self-adjoint. 
In particular, this implies the weak convergence of the corresponding spectral measures at the 
root and the compact convergence of their associated Cauchy-Stieltjes transforms : 



m Gn > m G and u Gn ' — > VG- 

n— voo n— yoo 

Note that this last statement does not depend anymore on the way G,G\, . . . have been embed- 
ded. We have thus established the following continuity result: 

Proposition 11 Let G, G%, G2, ■ ■ ■ be elements ofG* whose adjacency operators are self-adjoint. 
Let /xg,/UGd • • • denote the associated spectral measures at their root, and mc,mG 1 , ■ ■ ■ the cor- 
responding Cauchy-Stieltjes transforms. If G n > G, then 



m Gn > m G and ^i Gn > /x G . 

n— XX) n— yoo 

As a consequence, when G, G%, G2, ■ ■ ■ are random elements of G* , the same implication holds 
with all convergences being replaced by their distributional versions. More precisely, if the law 
of G n converges weakly to that of G, then 

V(H) , V{V(R)) 
mo n > m G and fx Gn > \x G . 



3.3 Connection with the empirical spectral measure of a finite graph 

In the case of a finite (non-rooted) graph G n = (V n , E n ) on n vertices, the adjacency operator A n 
is a particularly simple object: it is bounded and self-adjoint, and it has exactly n eigenvalues 
Xi(A n ) > . . . > X n (A n ) (with multiplicities), all of them being real. Moreover, £ 2 (V n ) = C™ 
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admits an orthonormal basis of eigenvectors (b\, . . . ,b n ), a priori different from the canonical 
orthonormal basis (e„) ve y n , such that 

n 

VxeC n , A n x = J2 x i( A n){x,bi)bi. 
1=1 

If (G n ,v) denotes the graph G n when rooted at v, the spectral measure at the root is simply 

n 

HG n ,v) = ^2\{ b ue v )\ 2 5 x . {An) . 
i=i 

In fact fi(G n ,v) can De interpreted as the local contribution of vertex v to the empirical spectral 
measure fi n of G n . Indeed, the above formula implies 

1 1 n 

~ Yl HGn,v) =~Y1 S *i(An) = Vn- (16) 

V&fn 1=1 

Note that the left-hand side can be reinterpreted as the expectation of [1(g„,0) under a uniform 
choice of the root 0. More generally, if G n is a random graph on n vertices, we denote by 
U(G n ) the random element of Q* obtained by rooting G n at a uniformly chosen vertex, inde- 
pendently of the random edge-structure. Similarly, we define [^(Gn) as the random element 
((G n ,0i), (G n ,02)) in Q* x Q*, where (0i,02) is a uniformly chosen pair of vertices. Finally 
we let [i n denote the (random) empirical spectral measure of the adjacency matrix of G ri . With 
this notation, 

Corollary 12 // U(G n ) converges weakly to a rooted GWT T whose degree distribution has 
a finite second moment, then 

lim Kfi n = EfiT, 

ri^>oo 

where [xt denotes the local spectral measure at the root of T. If moreover U2{G n ) converges 
weakly to (Ti,T 2 ), two independent copies of T , then in probability, 

lim \i n = E/iy. 

n— »oo 

In the above mentioned cases of Erdos-Renyi random graphs and random graphs with asymp- 
totic degree distribution F*, the assumption on U2{G n ) is easily checked. This corollary implies 
that the study of the limiting spectral measure of random tree-like graphs boils down to the 
study of the local spectral measure at the root of the limiting GWT. As we have seen, the latter 
is fully characterized by a simple RDE involving its Cauchy-Stieltjes transform. Note however 
that this result does not give the full statement of Theorem [TJi) : the almost sure convergence 
will be considered later. 

Proof of Corollary [TM By (I16p . we may write for any bounded continuous function /: M — > R, 
E / fdn n = ~J2 E f f d HGn,0) > E / f d VT, 
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where the convergence follows from the weak convergence U(G n ) — > T and the continuity result 
stated in Proposition [TTJ This is exactly saying that E/i n — > K^t- If moreover, U2(G n ) converges 
weakly to {T\,T2), then by the same argument, 

E ( / fdfl n ) = E ( l fMGn,0l) [ f d HGn,0 2 )) ~ ~> ( E f f^) , 

\Jr J n 0i£Vnj02€Vn \Jr Jr J V it ) 

and therefore, the second moment method sufficies to conclude that: 

/ fdfin P > E / fd/jlT, 
Jr n ^°° Jr 

which is exactly saying that /x n — > fix in probability. □ 



3.4 Main result : convergence of the rank 



We are now in position to state the main result of this paper. We consider a sequence of finite 
random graphs G\ , G2, ■ ■ ■ converging in distribution (once uniformly rooted) to a GWT whose 
degree distribution F* has a finite second moment. As above, 99* (x) = J2k F*(k)x k denotes the 
generating function of , and we consider the function 

M : x £ [0, 1] (-)• ip'^(l)xx + </j*(1 — x) + <£>*(1 — x) — 1, where x = — x)/^(l). 

Recall that M'{x) = (p"(l — x)(x — x) so that M(x) is a local extremum if and only if x = x. 

Theorem 13 Assume that U2(G n ) converges weakly to (T\,T2), two independent copies of a 
GWT whose degree distribution has a finite second moment. If the first local extremum of M 
is the global maximum, then in probability, 

lim — rank(yl n ) = 1 — max M(x). 

n->oo n x£[0,l] 

Moreover, a simple sufficient condition for the assumption on M to hold is that <p" is log-concave. 

If the assumption U2(G n ) — )• (T\,T2) is replaced by the weaker U(G n ) — > T, then we only have 
convergence of the expected rank. 

The log-concavity of ip" is a sufficient condition for the first local extremum of M to be a 
global maximum. Setting h : x 1— > x — x, we find 

u / , ,t„ , ^'(l-x)^'(l-x) , . , . , y'l{l-x)ip'l'{l-x) ip'l'a-x) 
Vx G (0, 1), h (x) = — t —7— — am, with g(x) = -7— — — — -. 

Now, if ip" is log-concave, then x \-> ip'"(x)/tp"(x) is non-increasing on (0, 1), and therefore, g is 
decreasing (as the difference of a decreasing function and a non-decreasing one) . Consequently, 
h" can vanish at most once on (0, 1), hence h! admits at most two zeros on [0, 1], and h at most 
three. The unique root x c of x = x is always one of them, and if x is another one, then so is x. 
Therefore, only two cases are possible : 
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- either x c is the only zero of h ; then h(0) > and h(l) < 0, so M is maximum at x c . 

- or h admits exactly three zeros x_ < x c < x + ; in this case the decreasing function g has 
to vanish somewhere in (0, 1), so h" is positive and then negative on (0, 1). Consequently, 
h is decreasing, then increasing, and then decreasing again. In other words, M is minimum 
at x c and maximum at 

In both cases, the first local extremum of M is its global maximum. 

The remaining part of this section is devoted to the proof of Theorem 1131 First, recall that 
n rank(.A n ) = 1 — /z n ({0}). From Corollary 1121 we have in probability, 

limsup^ n ({0}) <E^ T ({0}). 

n 

In order to prove Theorem 1131 it is thus sufficient to establish that 

liminfE / u ri ({0}) > max M(x). (17) 

n x6[0,l] 

To do so, we will use the Karp-Sipser leaf removal algorithm, which was introduced in [15] to 
efficiently build a matching (ie a subset of pairwise disjoint edges) on a finite graph. 

For our purposes, the leaf removal algorithm on a locally finite graph G = (V, E) can be 
described as an iterative procedure that constructs two non-decreasing sequences (At) t>0 and 
(Bi) t>0 of subsets of V as follows : we start with 

A = {v G V : deg{v; G) = 0} and B = 0. 

Thus, Aq is simply the set of isolated vertices in G. Then, at each step t G N, we let Gt be the 
subgraph of G spanned by the vertex-set V t = V\{At UBt U Pt), where Po = 0. We denote by 

L t = {v£V t :deg(v;G t ) = l} 

the set of its leaves. We also introduce the set of vertices that are adjacent to those leaves : 

W t = {ve V t \Lt : 3u G L t , (uv) G E}. 

We add to Pt the set of pairs of adjacent vertices in Lj: 

Pt+i = P t U{ve L t , 3u G L t , (uv) G E}. 

Then we set 

A t+ i = At U {u G L t : 3v G W t , (uv) G E} and B t+ \ = B t U W t . 

In words, for any leaf u of Gt whose (unique) neighbor v is not a leaf, we add u to At and v 
to Bt- Then to obtain Gt+\, all nodes in At+i U Bt+i U Pt+i ar e removed from G (note that to 
obtain Gt+i, all leaves from Gt are removed with their adjacent vertices). If Lt becomes empty, 
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we have (At+i, Bt+i) = (At,Bt), and the algorithm stops. Finally, in the case where the graph 
G is finite, we define 

LR t (G) = \A t (G)\ - \B t (G)\. (18) 

Note that for any finite graph G, the sequence (LKt(G))t>o is non-decreasing. Note also, that 
the leaf removal algorithm is well defined for a (possibly infinite) locally finite graph but the 
definition (|18|) makes sense only for finite graphs. The Lemma below states a connection between 
these numbers and the rank of the adjacency matrix of G. It was first observed in [UJ, and a 
proof can be found in [TJ] . 

Although we will not need it here, let us make for completeness the following observation, 
which was the original reason why this algorithm was introduced for finite graphs: each time 
a vertex v is added to Bt, one may arbitrarily associate it with one of its neighboring leaves 
u v £ At- Similarly, for every vertex v added to Pt, define u v as its other neighboring leaf in P t . 
The edge-set {(vu v ),v G Bt U Pt} is then a matching of G, and it is contained in at least one 
maximum matching of G. Since the graph is finite, the algorithm stops at a finite time t*. The 
subgraph of G spanned by the vertex-set V\(Af U Bt* U Pt*) is a graph with minimal degree at 
least 2 called the core of the graph. 

Lemma 14 For any finite graph G with adjacency matrix A, and any t G N, 

dimker(A) > |LR t (G)|. 

Proof. Let u% G Lq(G) be a leaf of G and v its unique neighboring vertex. Let G' = G\{u\, v} 
and A(G') the adjacency matrix of G' , we have 

dim ker .4(G) = dimkerA(G'), 

(see [6j). Now, if {u±, ■ ■ ■ , u a } C Lq{G), is the set of leaves adjacent to v, then {112, ••• , u a } are 
isolated vertices in G'. The vectors e U2 , ■ ■ ■ ,e Ua are thus eigenvectors of the kernel of A'. By 
orthogonal decomposition, we deduce that 

dimker^4(G) = a — 1 + dimker(A(G\{v,ui, • • • ,u a })). 

By linearity, we obtain that for any integer t, 

dimkerA(G) = \A t (G)\ - \B t (G)\ + dimker(^(G\(A U B t U P t )) > \A t (G)\ - \B t (G)\. 

□ 

The lower bound (]17p will now follow from the following proposition. 

Proposition 15 Let T be a rooted GWT whose degree distribution F* has a finite mean. Then, 

lim P(0 G A t {T)) - P(0 G B t (T)) = M(x ). 

t— >oo 

where xq G [0, 1] is the location of the first local extremum of M. 
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Proof of Proposition \lb\ The argument is close to that appearing in |15l Section 4] . For any 
vertex i 7^ 0, we run the leaf removal algorithm on Tj which is the tree T\ with an additional 
infinite path starting from i. We first compute the corresponding probabilities at = P(i G At{T\)) 
and fit = P(i G Bt{T\)). For our purpose, adding the infinite path amounts to increase artificially 
the degree of the root by 1: to be a leaf in Tj, the root needs to be isolated in T\. By construction, 
i is in Bt{T\) if and only if one of its children k is in At(Ty). Hence if N denotes the number of 
children of i, we have 

fit = E[l - (1 - a t ) N ] = l-(p(l- a t ), 

where ip is the generating function of N with distribution F given by (|1U|) . Similarly, i is in 
At{T\) if and only if all its children k are in Bt—^Ty), so that at = </?(A-i)- Hence for all t > 1, 
we have at = (p(l — (p(l — at-i)) and qo = 0. Since x i-> (p(l — tp(l — x)) is non-decreasing, at 
converges to a, the smallest fixed point of the equation x = ip(l — (p(l — x)), and fit converges 
to P = 1 — 99(1 — a). Note that ip(x) = </?*(x) / , where 93* is the generating function of F*. 
Hence, with the notation of Section \2A\ we have fi = 1 — a, a = a. In particular, we get xo = a. 

We now comupte P(0 G A t {T)) - P(0 G B t {T)). Recall that D(0) is the set of neighbors of 
the root 0. Here are all the possible cases : 

- if Vi G D(0), i G B t -x{f{), then G A t {T); 

- if there exists j G D(0)\(B t -i{f j )uA t {f i )) andVi G D(0)\j, i G B t -i(fi), then G A t (T); 

- if there exists i ^ j G D(0) such that i G and j £ £ t (Tj), then G 6 t (T). 

In all other cases, ^ •At(T') U Bt(T). In summary, we have 
P(0GA(T)) = p(viG J D(0),iGS t _ 1 (f i ) N 



+ P (3j G D(0)\(S t _i(7}) U A(Tj)), Vi G £>(0)\j, i G Bi-i(fi)) 
^(A-i) + (l-A-i-aty.(A-i), 



»(0GB t (T)) = P(3i^jGL>(0), i€^(r J ),j^B t _i(r J ) 

3i G D(0), i G A(Ti)) - P (3i G D(0), i G A(T ; ), Vj G D(0)\i, j G ^-i(Tj) 



Hence, 



1 - - o t ) - a t tpi(Pt-i)- 



limP(0G At(T))-F(0eB t (T)) = </>*(/3) + (1 - P)cp'M + <p.(l - a) - 1 

t— >oo 



M(a) = M(x ), 



where we have used the identities : /3 = 1 — a, (p'*(x) / (p'*(l) = 1 — x and 1 — /3 = a. □ 
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Proof of Theorem 1 13\ As already pointed out, it is sufficient to prove (|17p . From Lemma [HI 
for any integer t, 

E/i„({0}) > -ELR t (G n ) = P(0 G A(G„)) - P(0 € B*(G n )), 
n 

where is the uniformly drawn root of U(G n ). Note that the events {0 G .At(G n )} and 
{0 G ^(G n )} belong to the a— field generated by (G n ,0) t +i- Thus the convergence of U(G n ) 
implies that for any t G N, 

lim P(0 G A(G n )) - P(0 G B t (G n )) = P(0 G A(T)) - P(0 G B t (T)), 

n— >oo 

where T is a rooted GWT with degree distribution F* (this is a standard application of the 
objective method |3j) . We then use Proposition [15] to conclude the proof of Theorem [T3l □ 



4 Conclusion 

As explained in the Introduction, the condition on M in Theorem [13] is restrictive and the 
convergence of the rank when this condition is not met (as in the example described in the 
Introduction) is left open. Without any condition on the function M, our work gives only the 
following bounds: Assume that LT 2 (G n ) converges weakly to (Ti,T2), two independent copies of 
a GWT whose degree distribution has a finite second moment, then in probability, 

1 1 

1 - max Mix) < liminf -rank(^ n ) < lim sup -rank(A n ) < 1 - M(x Q ), (19) 

X6[0,l] n->oo n n^oo U 

where xq is the first local extremum of M. For example, if the sequence of graphs converges 
weakly to a GWT with degree distribution F* with F*(l) = 0, i.e. with no leaf, then xq = and 
M(0) = -F*(0) so that the upper bound in (|19j) is trivial. 

Our proof for the upper bound on the rank of A n relies on the analysis of the leaf removal 
algorithm on the graph G n . As explained above, this algorithm when applied to a finite graph 
produces a matching and a subgraph of minimal degree 2 called the core. It turns out that 
the RDEs (I12p and (|13j) also appear in the the analysis of the size of maximal matchings on 
graphs [TTJ. In particular, if the size of the core is o(n), the leaf removal produces an (almost) 
maximal matching (with error o(n)) and the bounds in (|19f) match. If the size of the core is not 
negligible but the bounds in (|19p match (as for example in the case where ip" is log-concave), 
our result shows that the asymptotic size of the kernel of the core is zero. In [11] , it is shown 
that this case corresponds to the situation where there is an (almost) perfect matching on the 
core of the graph. However as soon as M(xq) ^ max^mi] Mix), for any maximal matching, 
there is a positive fraction of vertices in the core that are not matched [UJ . In this latter case, 
the convergence of the rank is left open. 
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5 Appendix: proof of Theorem [I] 

In the case where is the Poisson(c) distribution, we simply have 

Vx G (0, l),tp(x) = <p*(x) = exp(c(x — 1)), 

whose second derivative is clearly log-concave. We may therefore apply Theorem [13] to the 
sequence of Erdos-Renyi graphs (Gn)neN- 




Figure 2: From left to right : plot of M for c = 2, c = e and c = 3. 

To complete the proof of Theorem[Tl it only remains to improve the convergence in probability 
into an almost sure convergence. This is performed by a standard exploration procedure of the 
edges E n of the graph G n . For 1 < k < n, we define the random variable in {0, l} k : 

Xk = (Aik)i<i<k- 

By construction, the variables (Xk)i<k<n are independent random variables. Note also that 
the upper half of the adjacency matrix A n is precisely {X\, ■ ■ ■ ,X n ) and we may safely write 
A n = A(X 1 ,--- ,X n ). 

For 1 < i < n, let Ai(X\, ■ ■ ■ , X n ) be the principal minor of A obtained by removing i-th row 
and column. If Ai < • • • < X n and A^j < • • • < \ n -i,i denote the eigenvalues of A(X±, ■ ■ ■ ,X n ) 
and Ai{X\, ■ ■ ■ , X n ), by the Cauchy interlacing theorem, for all 1 < j < n — 1, 

Aj < Xj.i < Aj-fi. 
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In particular, 

|dimkerA(Xi, ■ ■ ■ , X n ) - dimkevA i (X 1 , ■ ■ ■ ,X n )\ < 1. 

We note that • • • , X n ) does not depend on X,- t . Therefore, for all (xj € {0, l}- 7 ), 1 < j < n, 

^€{0,1}*: 

|dimker^(xi, • • • , Xi-i, x iy x i+ i, ■■■ ,x n ) - dimkerA(xi, • • • , Xi-i, x\, x i+ i, ■ ■ ■ ,x n )\ < 2. 

In other words, the function (xi, ■ ■ ■ , x n ) (->• dimker^4(xi, • • • , x n ) is 2-Lipschitz for the Hamming 
distance. By a standard use of Azuma's martingale difference inequality we get 

( -t 2 ^ 

F(|dimkerA(Xi,-- - ,X n ) - Edimker,4(Xi, ■ ■ ■ ,X n )\ >t)< 2exp 

\ 8n 

^From Borel Cantelli Lemma, we obtain that almost surely, 

dim ker A(X X , • • • , X n ) - Edim ker A(X\ , • • • , X n ) 
lim = 0. 

n n 

Since we have already proved that Edim ker -A(Xi, • • • ,X n )/n converges to max^^i] M(x), we 
deduce that dimker^4(Xi, • • • ,X n )/n convergences a.s. to max^Q ^ M(x). 

It remains to deal with the almost sure convergence in Theorem [JJ(i). We have already proved 
that [i n converges in probability to [i. Henceforth E/x n converges to [i. It is thus sufficient to 
prove that almost surely, for all t £ R, /x n ((— oo,i]) — E// n ((— oo, t]) converges to 0. The next 
lemma is a consequence of Lidskii's inequality. For a proof see Theorem 11.42 in [4] 



Lemma 16 (Rank difference inequality) Let A, B be two n x n Hermitian matrices with 

i(A) and hb = ^ELi^OB)- 



empirical spectral measures ha. = h Ya=i $\(A) an d Hb = ^ ^(B)- Then 



sup \ha({— oo, t]) — oo, t])\ < — rank(^4 — B). 

tm n 

Again, we view fj, n as a function of (X x , ■ ■ ■ ,X n ), and write ii n = Hrx u — ,x n )- Note that for all 
(xj G {0, 1} J '),1 <j < n, x\ € {0, 1}\ A(xi, ■ ■ • , x;_i, Xi, x i+1 , ■ ■ ■ , x n )-A(x l , ■ ■ ■ , x\, x i+1 , ■ 
has only the i-th row possibly different from 0, and we get 

rank (A(x\, • • • , Xi—i, x^, ■ ■ ■ , x n ) — A(x\, ■ ■ ■ , Xi—i, x^, • • • , £ ra )) < 2. 

Therefore from Lemma [TBI for any real t, 

1 



< 
n 



M(a5l,— ,SDi-l,Xi,Xi+l,— ,X n ) (( °°, t]) fXi-ltX'iJCi+l,— ,X n ) (( °°>^]) 

Again, Azuma's martingale difference inequality leads to 

— 

We deduce similarly from Borel-Cantelli lemma that fi n ((—oo,t]) — E/%((— oo, t]) converges a.s. 
to and the proof of Theorem [1] is complete. 
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